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Abstract 

We calculate leading-order medium photon yields from a quark-gluon plasma using (3-1-1)- 
dimensional anisotropic hydrodynamics. Non-equilibrium corrections to the photon rate are taken 
into account using a self-consistent modification of the particle distribution functions and the corre¬ 
sponding anisotropic hard-loop fermionic self-energies. We present predictions for the high-energy 
photon spectrum and photon elliptic flow as a function of transverse momentum, shear viscosity, 
and initial momentum-space anisotropy. Our findings indicate that high-energy photon produc¬ 
tion is sensitive to the assumed level of initial momentum-space anisotropy of the quark-gluon 
plasma. As a result, it may be possible to experimentally constrain the early-time momentum- 
space anisotropy of the quark-gluon plasma generated in relativistic heavy-ion collisions using 
high-energy photon yields. 
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I. INTRODUCTION 


The goal of the ongoing relativistic heavy-ion collision experiments at Brookhaven Na¬ 
tional Laboratory’s Relativistic Heavy Ion Collider (RHIC) and at CERN’s Large Hadron 
Collider (LHC) is to produce and study the properties of the quark-gluon plasma (QGP). It 
is now well accepted that a few microseconds after the Big Bang the entire universe consisted 
of an extremely hot and dense QGP. In order to reproduce these conditions terrestrially, rel¬ 
ativistic heavy-ion collisions are used. As a result, extremely high temperatures and energy 
densities are created within a very small volume (~ 4000 fm^). At these high temperatures, 
quarks and gluons no longer remain conhned within nucleons and one instead generates 
deconhned nuclear matter called a QGP. Immediately after the initial nuclear impact, the 
QGP generated in relativistic heavy-ion collisions cools by expansion and, below a certain 
(pseudo-)critical temperature (T^ ~ 165 MeV), the quarks and gluons recombine to form 
hadrons. After the transition to hadrons, the system may undergo further expansion and 
cooling before full chemical and kinetic freeze-out takes place. The resulting particle produc¬ 
tion and associated radiation from the event are then analyzed in order to infer information 
about the properties of the QGP. 

One of the key outstanding questions in the study of the QGP is the question of the time 
scale for the thermalization and isotropization of the matter created in relativistic heavy-ion 
collisions. Theoretical calculations in both the weak-coupling and strong-coupling limits hnd 
that the QGP created immediately after the initial nuclear impact (r ~ 0.2 fm/c) is highly 
anisotropic in local rest frame (LRF) momentum, however, there are currently no clear 
experimental observables that can be used to conhrm this expectation and constrain the 
degree of early-time momentum-space anisotropy. For this purpose, radiation of photons 
and dileptons are promising signals since they can be used to probe the initial state of 
heavy-ion collisions. Unlike hadrons, which are emitted from the freeze-out surface after 
undergoing intense re-scatterings, photons emerge from all phases of the expanding hreball; 
initial hard scatterings, pre-equilibrium phase, near-equilibrium phase, and hadronic phase. 
Since photons are electromagnetic probes, they interact only weakly with the QGP (a a^) 
and their mean free path is much larger than the typical system size (~ 10 fm). As a result, 
once produced, they do not undergo signihcant interactions with the medium and carry 
largely undistorted information about the circumstances of their production to the detector. 
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For the most part, in the past calculations of the QGP photon production rate have 
been performed assuming a perfectly thermalized, weakly-coupled QGP when using hydro¬ 
dynamics for the background evolution. Within this framework a complete calculation of 
the thermal photon rate at 0{e‘^gs‘^) has been available for a decade [1] and the next-to- 
leading-order (NLO) correction 0{e^gs^) to thermal rate has been computed recently [2]. 
At low-temperatures, below the pseudo-critical temperature for the QGD phase transition, 
where dense QGD matter can be modeled as a hadron resonance gas, effective Lagrangian 
approaches have been used, see e.g. Ref. [3]. The success of viscous hydrodynamics applied 
to heavy-ion collisions [4, 5] suggests that it might be reasonable to assume that the medium 
is close to being in local thermal equilibrium. However, nonzero values of the QGD transport 
coefficients, resulting from nonzero mean free paths of the constituents, lead to deviations 
from local thermal equilibrium which increase with the local expansion rate. 

For example, in a dissipative QGP, a hnite shear viscosity causes the momentum distri¬ 
bution in the LRF to become highly anisotropic at early times, with the distribution falling 
off more rapidly in the directions in which the system expands (longitudinal cooling). In 
the past, various attempts have been made to determine the effect of viscous corrections 
on the photon emission rates in a QGP [6-8]. However, these previous works have a key 
shortcoming: They include only the viscous corrections to the local momentum distribution 
functions for the incoming and outgoing particles, but ignore viscous medium modihcation 
of the collision matrix element itself. For scattering processes in which the inclusion of 
medium effects is essential, one must self-consistently include viscous corrections. For ex¬ 
ample, when dynamical mass generation for the medium constituents serves as a regulator 
for IR divergences, viscous corrections to the distribution functions can lead to signihcant 
modihcations of the screening mechanism and, therefore, to the collision matrix element it¬ 
self. This problem was hrst addressed using a spheroidal form for the LRF momentum-space 
anisotropy in Ref. [9] with details and followup studies presented in Refs. [10-13]. 

Recently it has been shown how to extend the methods used in Ref. [9] to include the full 
shear-stress tensor modihcation of the one-particle distribution function [14-17]. The calcu¬ 
lation was based on photon production including Gompton scattering and quark-antiquark 
annihilation at leading order in ag. Extending the proof presented originally in Refs. [9, 10], 
Shen, Paquet, Heinz, and Gale were able to show that the Kubo-Martin-Schwinger relation 
holds in the hard-loop (HL) regime for any particle momentum distribution function that is 
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reflection symmetric.^ Using this, they were able to compute the hard and soft contributions 
to the rate separately, taking into account the modifications to the (anti-)quark self-energy 
necessary to make it consistent with a standard Grad-14 6f modification to the one-particle 
distribution function. 

One potential problem with using the standard Grad-14 form for the viscous correction 
to the particle distribution function is that, in the integrals that determine the photon 
production rate, one is integrating over all momenta. At high-momentum, the standard 
viscous approximation to the particle distribution function is not reliable and, as a result, one 
finds regions in the integration domain where the viscosity-corrected distribution function is 
negative. In this paper, we use instead a spheroidally-deformed form of the non-equilibrium 
distribution function that is positive by construction. We concentrate on photon production 
from the deconfined phase of the QGP’s lifetime and extend prior results performed in 
Refs. [10, 11, 18] to include a more realistic QGP background evolution. In order to improve 
upon previous works [10, 11, 18], in which a simple analytic model with an adjustable 
isotropization time was used for the QGP evolution, in this paper we use leading-order 
(3-1-1)-dimensional anisotropic hydrodynamics ((3-|-l)D aHydro) [19-22] to describe the non¬ 
equilibrium QGP evolution.^ Herein, we use the aHydro equations obtained from the zeroth 
and first moments of the Boltzmann equation with the collisional kernel treated in the 
relaxation-time approximation. The resulting dynamical equations describe the full (3-|-l)D 
spatiotemporal evolution of the transverse momentum scale A and spheroidal momentum- 
space anisotropy parameter ^ [24, 25]. The (3-|-l)D framework allows A, and the associated 
flow velocities to depend arbitrarily on the transverse coordinates, spatial rapidity, and 
longitudinal proper-time, however, herein we restrict ourselves to smooth Glauber-like initial 
conditions. 

We also mention that there there have been recent studies that have used various types of 
partonic and hadronic kinetic transport codes to address the problem of photon production in 
heavy-ion collisions [26-29] . Such codes automatically take into account the non-equilibrium 
quark and gluon phase space distributions, but they ignore the effect of non-equilibrium 
(anisotropic) screening in the problem. Typically in the kinetic approaches it is assumed 
that there is a cutoff on u- and t-channel exchanges which is either fixed or dynamically set 

^ One must also require that fq = fq, which was implicit in their proof. 

^ For a recent review of anisotropic hydrodynamics see Ref. [23] . 
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by the local Debye mass determined from the parton density. As shown in Refs. [9, 10], in 
order to properly regulate the infrared divergences in the photon production calculation, it 
is necessary to revisit the calculation of the of screening and use a quark self-energy which 
is consistent with the non-equilibrium (locally anisotropic) nature of the QGP. That being 
said, it may be that, within the accuracy required for QGP phenomenology, it is sufficient to 
simply have an isotropic infrared cutoff that does not take into account the non-equilibrium 
nature of the quark distribution function. In such a case, the transport codes would provide 
a quite reasonable calculation of photon production. 

The study presented in this paper is similar in spirit to other studies of photon pro¬ 
duction using second-order viscous hydrodynamics. Our goal is not to produce the most 
complete calculation of photon production from all possible sources and including fluctu¬ 
ating initial conditions etc., but to instead study the effect of self-consistently including 
the non-equilibrium modihcations of the quark distribution function and to systematically 
investigate the dependence of the resulting spectra and elliptic flow on the assumed shear 
viscosity and initial momentum-space anisotropy. Our work goes beyond prior viscous hy¬ 
drodynamics studies by linearizing around anisotropic background and, as a result, we are 
able to better describe early-time photon production and photon production near the trans¬ 
verse and longitudinal edges of the QGP. In addition, high-momentum photon production 
is treated in a more reliable manner since the anisotropic one-particle distribution function 
used to compute the photon rates is positive dehnite at all points in momentum space. Our 
results indicate that high-energy photon production is quite sensitive to the assumed level 
of initial momentum-space anisotropy of the QGP. As a result, it may be possible to ex¬ 
perimentally constrain the early-time momentum space anisotropy of the QGP generated in 
relativistic heavy-ion collisions using the high-energy photon spectrum. 

The structure of the paper is as follows. In Sec. II, we discuss the calculation of the 
photon rate from an anisotropic QGP. In Sec. Ill, we discuss how the rate is integrated over 
space-time to obtain the hnal differential photon yields. In Sec. IV, we specify the details 
of the hydrodynamical evolution utilized. In Sec. V, we present our numerical results for 
different initial conditions and different values of the shear viscosity to entropy density ratio. 
In Sec. VI, we present a discussion of our results, conclusions, and an outlook for the future. 
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II. PHOTON RATE IN ANISOTROPIC PLASMA 


In this section, we review the calculation of the photon production rate in an anisotropic 
QGP. This method was introduced originally in Ref. [10], but we review it here for complete¬ 
ness. To proceed, one separates the contributions to the rate into those corresponding to 
hard-momentum and soft-momentum exchanges which we detail separately below. In both 
cases, we take the quarks to be massless since, in the high-temperature and high-energy lim¬ 
its, the masses result in very small corrections to the relevant cross sections. In this paper 
we assume that, when written in terms of LRF momentum, the one-particle distribution 
function is of spheroidal (Romatschke-Strickland) form [24, 25], 


/.(k,{, A) = /=‘>(yk2+{(k.n)VA), 


( 1 ) 


where i = {g, g, g}, A is the transverse momentum scale, h is the direction of the anisotropy, 
and — 1 is a parameter reflecting the strength and type of the anisotropy. Above ^ and 

A should be understood to be helds that depend on both space and time. The function 
/eq is an equilibrium distribution function. In the following, we will suppress the explicit 
dependence of the anisotropic distribution function on ^ and A. 

1. Hard Contribution 

For the hard contributions, one can simply compute the Feynman diagrams corresponding 
to the Compton and annihilation processes. The rate of photon production due to in-medium 
quark annihilation can be expressed as 



?/ f 

x5\K^ + K2-Q-K:i) - + - 

t u 


( 2 ) 


where the Mandelstam variables are dehned hy t = {Ki — QY and u = {K 2 — Q)^, Cq^ = 2/3,^ 
f/^ = f d^k/(27r)^, and fq^g are the in-medium quark and gluon distribution functions. Note 
that herein capital letters, e.g. iFi, indicate four vectors. Henceforth, we also assume that 


^ Herein we include contributions from up, down, and strange quarks since these dominate the bulk of the 


QGP. 
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the system is charge-conjugation symmetric such that the distribution functions for quarks 
and anti-quarks are the same, i.e. fg = fg. 

In order to regulate the IR divergence associated with this graph, we hrst change variables 
in the hrst integration to P = Ki — Q and introduce an IR cutoff p* on the integration over 
the exchanged three-momentum p [30]. Here we choose spherical coordinates with the 
anisotropy vector h dehning the 2 ;-axis and we exploit the azimuthal symmetry around the 
z-axis to choose q to he in the x-z plane. Using the energy-momentum conserving delta 
function and expanding out the phase-space integrals gives 

2 /» rvn /»1 /»97r n / . \ /»r>o /*1 

dkk d{cos 6 k) 




d^q 


2 tt^ 


W [ °°dpp‘ f d(cos9,) rd4,,tdtP±^ 

~^Jp* J-1 Jo |p + q| 

X/q(k)[l + fg{p + k)]x-^/20(x) 


, (3) 


— 4^i 


with t = oj'^ — p‘^, u = {k — qY — (k — q)^, and uj = |p + q| — O'. The azimuthal angles <pi are 
dehned through 

uj‘^ — p'^ + 2 k{pj — p cos 6 p cos 6 k) 


cos(0i - 0p) = 


2 pk sin 6 p sin 6 k 


(4) 


and X > 0 is given by 

X = Ap^k"^ sin^ 6 k sin^ 6 p — \uY — p^ + 2k{uj — p cos 6 p cos 6 k)] ^ . (5) 


The rate of photon production from the Compton scattering diagrams can be obtained from 


pdl^ = -USTT^elasa 


d^q 


f /.(ki) 

f /.(k2) f 

1 - /.(ks) 

/ki 

ik2 ^2 A3 

h 


x5\K, + K 2 - Q - K 3 ) 


s t 

Us 


( 6 ) 


where the Mandelstam variables are dehned by s = {Ki -|- ^" 2 )^ and t = {Ki — Q)^. After 
changing variables to P = Ki — Q and continuing as with the annihilation process one 
obtains 

2 /»r5n /»! /»27r 


^dRcom OtgOl. 


dP'q 


2 tt^ 


^ poo pi 

/ dpp"^ / 
i=i -^P* -^-1 


d cos 6 . 


d(f). 


fqip + q) 


^00 




p I u.gyp- -— f dkk dcos6k 

0 |P + q| Jo J-i 


x/g(k)[l-/g(p + k)]x ^''^0(x) 


s t 


, (7) 


4^k — 4^i 
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with t = and s = (cn + g + fc)^ — (p + k + q)^. The total photon prodnction rate 

from hard processes is given by the snm of Eqs. (3) and (7). 



(Pq 



( 8 ) 


We nse Vegas Monte-Carlo integration [31] to evalnate the remaining hve-dimensional in¬ 


tegrals in Eqs. (3) and (7). The total hard contribntion (8) has a logarithmic infrared 


(IR) divergence as p* —)■ 0. This logarithmic IR divergence is cancelled by a corresponding 
nltraviolet (UV) divergence in the soft contribntion which we will describe next. 

2. Soft Contribution 

Next we calcnlate the contribntion involving soft-momentnm exchange. We refer the 
reader to Ref. [9] for farther details of the calcnlation. We nse the Keldysh formnlation, 
which is appropriate for non-eqnilibrinm systems [32]. The components (12) and (21) of the 
polarization tensor are related to the emission and absorption probability of the particle 
species nnder consideration, respectively [33-35]. Dne to the very low rate for photon 
absorbing back reactions, the rate of photon emission can be expressed as [36] 



(9) 


from the trace of the (12) element of the photon polarization tensor. We evalnate ( 1112 )^^ 


nsing the HL resnmmed fermion propagator derived in Ref. [9] wherein the anthors demon¬ 


strated that the needed off-diagonal components of the fermion self-energy can be expressed 
in terms of the retarded self-energy 



( 10 ) 


where 


/(k)^2(/,(k) + /,-(k))+4/,(k). 
Taking the HL limit where appropriate, one hnds that 


( 11 ) 



( 12 ) 


where 










- Total (Numerical) 

- Soft (Numerical) 

- Hard (Numerical) 

- Total (Analytic) 

. Soft (Analytic) 

- Hard (Analytic) 


FIG. 1. Dependence of the isotropic photon rate on the separation scale p*. The various curves 
show the total rate determined numerically, the soft and hard rates determined numerically, the 
LO rate analytic result of Braaten and Yuan [30], and the separate Braaten and Yuan results for 
the soft and hard rates. 

The tensor A is defined throngh 

A = rpil-h^YTL rui 

af-A ) (p-s(p))^ ^ Aa-S"(P))2’ ' * 

where a star indicates complex conjngation. To compnte the soft photon rate (9), we nn- 
merically evalnate Eqs. (10) and (12) with an UV cntoff p* placed on the length of the 
three-momentnm in Eq. (12). 

In Fig. 1 we show the soft and hard contribntions to the isotropic (^ = 0) photon rate and 
the total rate obtained by snmming these contribntions. We compare onr nnmerical resnlts 
with the analytic estimates of Braaten and Yuan [30] for both the individual contributions 
and the total rate. For this figure, we consider the case of a realistic coupling with as = 0.32 
and a photon momentum which is hve times the corresponding transverse momentum scale 
A.^ As can be seen from Fig. 1, there is a logarithmic IR divergence at small p*/A and 
logarithmic UV divergence at large p*/A. However, even for large there is a window 
over which the result does not depend strongly on the choice of p*. In practice, we use the 
principle of minimal sensitivity (PMS) to set p* by requiring the derivative of the rate with 
respect to p* to vanish (minimum of the black solid curve in Fig. 1). The other things we 

Since ^ = 0, A can be identified with the temperature T in this figure, but we have kept the labels general 
for comparison with other figures. 
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see from this figure are that the analytic calculations do well in capturing the asymptotic 
regimes of each contribution, but that the total result obtained from analytic method is 
lower by approximately a factor of 2 due to the fact that the analytic approximations made 
to obtain the Braaten-Yuan rate fail to accurately reproduce the numerically integrated 
rate.® 

In Fig. 2 we plot the photon rate as a function of the scaled photon energy E/K for (a) 
^ G {0,10} for j/ = 0 and (b) .^ = 10 for ?/ G {0,0.88,20}. In both panels of Fig. 2, the 
lines are the result of evaluating the rate at the PMS value of p* and the shaded regions 
indicate the variation of the result obtained when varying p* —)■ 2p*. In the high-energy 
and weak-coupling limit, the dependence on p* formally vanishes, however, for realistic 
couplings, direct numerical evaluation of the integrals defining the rate allows us to gauge 
this uncertainty. From the results we hnd that there is a ^ 30% variation of the photon 
rate at E/K = 1 and a ^ 20% variation at E/K = 10. From Fig. 2(b) we conclude that 
there is a signihcant rapidity dependence of the photon rate when the system is momentum- 
space anisotropic, with production peaked at mid-rapidities. This is simply due to the 
dominance of forward scattering coupled to production from a momentum-space anisotropic 
source with {p^) > {p\). Note, however, that, if the system were exactly boost-invariant 
then, when integrated over all spatial rapidity, the final photon spectrum would not depend 
on the photon rapidity. If, on the other hand, the system is not boost-invariant, then we 
expect to see a suppression of photon production at forward/backward rapidities if the LRF 
distribution is oblate. 

III. PHOTON SPECTRUM 

As mentioned previously, in this paper we want to study the impact of space-time de¬ 
pendent anisotropies on the photon differential spectrum at high transverse momentum. In 
this way we can probe the early stages of the QGP, when the anisotropies are expected to 
be the largest. In order to accomplish this, we need to convolute Eqs. (8) and (9) with the 
space-time dependence of A and ^ to integrate the rate over the QGP space-time volume 
and obtain the hnal spectrum of photons emitted from the deconhned QGP over its lifetime. 
® If one takes very small Og, e.g. Og = 0.01, then the numerical and analytic results agree nearly perfectly 
for large g/A, giving us confidence in our numerical methods. 
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FIG. 2. Rate as a function of the scaled photon energy for (a) ^ G {0,10} for y = 0 and (b) ^ = 10 
for y G {0,0.88,20}. The shaded bands show the change in the result obtained by varying the 
separation scale p* —)■ 2 p*. 


We now discuss how this is implemented. 

The photon four-momentum is parametrized in the standard manner 

= q± (cosh y, cos (pq, sin (j)q, sinh y ), (15) 

where y = ln[(iil -|- q\\)/{E — q'||)]/2 is the momentum-space rapidity. Above g_L, q\\, and (j)q 
are the transverse momentum, longitudinal momentum, and momentum azimuthal angle, re¬ 
spectively. For space-time, we use the usual Milne hyperbolic parametrization for describing 
heavy-ion collisions within the relativistic hydrodynamics framework 


= (r cosh <^, x_i_, r sinh <^) . 


(16) 


In Eq. (16), r = \/t‘^ — is the longitudinal proper time, is a two-vector contain¬ 
ing the transverse coordinates, and ? = tanh“^( 2 ;/f) is the space-time rapidity. With 
these parametrizations, the differential measures for four-momentum and space-time are 
d^q = E dy q±dqj_d(j)q and d'^x = r dr dq d‘^x±, respectively. This allows us to calculate the 
differential spectrum using 


dN 


q±dq±dy 


I'2'k 

= / d(l)q 

Jo 


, „dR 


(17) 


LRF 


The integration over the space-time volume is performed solely in the deconhned QGP stage. 
We only include contributions from regions that have an effective temperature higher than 
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a critical temperature, i.e. Teg = > Tc with 7^(0 defined in Eq. (47). In all results 

shown herein, we assume Tc = 175 MeV. We will gauge the sensitivity of our results to this 
assumption in Sec. V. We will assume that when the system reaches Tc, all QGP medium 
emission stops. We do not consider the emission from the mixed/hadronic phase herein. 

Finally, we mention that we evaluate the left-hand-side of Eq. (17) in the center of mass 
of the colliding nuclei (LAB frame), while the photon rate is calculated in the LRF of the 
emitting region. Therefore, before evaluating Eq. (17), we have to boost the LAB frame 
momentum to the LRF of the fluid cell using a Lorentz boost ~ ^ where the 

Lorentz boost tensor 



( 7 


-l^y 

-IVz 

-7Vx 

l + (7-l)S 

(7-1)^ 

(7-1)"^ 

-IVy 

(7-1)^ 

l + (7-l)S 

(7-1)^ 


(7-1)^ 

(7-1)^ 

l + (7-l)S 


( 18 ) 


depends on the four-velocity of the fluid element u^{x^) = ■j{l,Vx,Vy,Vz), with 7 = 
1/Vl — and V = + Vy + v^. Making use of Eqs. ( 8 ) and (9) in Eq. (17), we obtain 

the photon spectrum including the effect of a space-time-dependent momentum anisotropy 
and taking into account the effect of the dynamically-generated collective flow of the QGP. 


IV. (3-hl)D ANISOTROPIC HYDRODYNAMICS 

As mentioned above, in order to obtain the predictions for the differential photon spec¬ 
trum expected to be produced from the QGP phase, we must integrate over the full 
space-time history of the QGP. For this purpose we use a (3-|-l)D leading-order spheroidal 
anisotropic hydrodynamics (aHydro) code [19, 20, 37-41]. The aHydro framework reduces 
to second-order viscous hydrodynamics in the limit of small anisotropy [ 21 ], but reproduces 
the dynamics of the QGP more reliably when there are large momentum-space anisotropies. 

We assume that the QGP created during the collision of the heavy ions evolves through a 
non-equilibrium state and that the quark and anti-quark one-particle distribution functions 
are well-approximated by a time-evolving distribution of the form specified in Eq. (1) both 
at early times and late times. At the same time, we also assume that, although the system 
is highly anisotropic, it may still be, to good approximation, described using hydrodynamic- 
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like degrees of freedom, such as energy density and pressures.® As a result, the detailed 
microscopic description of the system can be replaced by an effective description written in 
terms of simple physical laws, such as conservation of energy and momentum. 

The equations of motion for the anisotropic system are obtained by starting from kinetic 
theory, assuming that the form of the distribution function of the system is known at the 
leading order and given by the form (1). This can be performed by taking moments of 
the Boltzmann kinetic equation. As usual, the collisional kernel must be specihed in the 
Boltzmann kinetic equation. Here we will assume that the collisional kernel can be treated 
in the relaxation-time approximation (RTA) such that 

p"a,/ = -^(/-/eq), (19) 

where Teq is the microscopic relaxation time which can depend on position and time. Taking 
the hrst moment of the Boltzmann equation results in the energy-momentum conservation 
equation 

= 0 . ( 20 ) 

Taking the zeroth moment of the Boltzmann equation results in the particle production 
equation 

= -u/- -^ . (21) 

^eq 

At leading order, the aHydro energy-momentum tensor has the form typical for a spheroidally 
anisotropic system 

= (e + P^) uV - - (Px - P\\)z^^z^ , (22) 

and the particle flux is dehned in the standard manner^ 


f = nu^. (23) 

In Eqs. (22) and (23) e, n, T||, and Pi stand for energy density, particle density, longitudinal 
pressure, and transverse pressure, respectively. The four-vector is orthogonal to and in 
the LRF points in the longitudinal direction (identihed with the direction of the dynamically- 
evolving anisotropy in the system, h) [40]. 

® This assumption has been tested elsewhere by comparing the predictions of anisotropic hydrodynamics 

to exact solutions of the Boltzmann equation in a variety of special cases [42-49]. 

^ We assume vanishing chemical potential gradients. 
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Equations (20) and (21) provide a set of five independent partial differential equations 


DuS — — (e + Pi) 6u + (Pi — P\\) UyD^z^ , 
= (Pi - P||) + + Pi) , 


Pu^l. 


U± 

£ + Pi 


Ui ■ ViPi 



+ DuP± + (Pi 


P\\)u^D^z’' , 


D 


U 



1 

ul{e + PA_) 


{Uxdy 


Pi. •) 


(24) 

(25) 

(26) 
(27) 


and 


Dui 

2(1 + 0 


3P„A 

A 





(28) 


respectively, for five parameters: the three independent components of the four-velocity 
the transverse momentum scale A, and the anisotropy parameter In the above equations, 
we use boldfaced letters with a + subscript to indicate two-dimensional vectors in the trans¬ 
verse plane, e.g. ui = {ux^Uy) and Vi = {dx,dy). We have also introduced a compact 
notation for the convective derivative Du = u^dy,, the longitudinal derivative = z^dy, 
and the expansion scalars 9u = dyU^ and 9z = dyZ^. 

We use the following parametrizations of the LAB frame four-velocity of the fluid and 
the space-like four-vector z'^ 


= (uo cosh'd, ui, Mo sinh'd), 

(29) 

z^ = (sinh -d, 0, cosh'd), 

(30) 

where we introduced the longitudinal rapidity of the fluid cell '9. 

normalization condition, u^Uy = 1, one has 

Using the four-velocity 

Uq = sjl + u\, 


= ^Jul + Ul- 

(31) 

With the parametrizations (29) and (30), one may calculate the following quantities appear- 
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ing in Eqs. (24)-(28), 


Du = u±-V±+ uqLi , (32) 

Ou = 'V± ■ u_i_ -|- LiUq + Uf)L 2 'i!^ , (33) 

Dz = L 2 , (34) 

0z = L,^, (35) 

u^D^z"' = UoL2-d , (36) 

ZyDuu'" = -Mo (u± ■ V± -h uoLi^ 10 , (37) 

where the two linear differential operators, Li and L 2 , are given by 

'' d 

Li = cosh(<^ — i0)dr — sinh(<^ — {})—, (38) 

T 

^ Q 

—L2 = sinh(^ —'d)^,-— cosh(<^ —'d)—. (39) 

r 


We also use the relation between the relaxation time Xeq and the shear viscosity to entropy 
density ratio fj = rj/s [19], 

= (40) 


A. Anisotropic equation of state 


In this paper we consider a system that consists of massless particles described by the 
anisotropic distribution function (1). Using standard kinetic theory definitions 


= 




(41) 




(42) 


where dK = d^k/ [(27r)^fc°], and the tensor decompositions specified in Eqs. (22) and (23), 
one can calculate the thermodynamic properties of the system [50] 


’^(A,0 


^iso(A) 

TttI’ 


£(A,0=77(e)Uso(A), 


Px(A,0=77x(Oao(A), 


P|l(A,0=77||(0Piso(A), 


(43) 

(44) 

(45) 

(46) 
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where riiso, ^iso, and Piso are the isotropic particle density, energy density, and pressnre, 
respectively, and 


no 

KiK) 

KiiK) 


1 

1 ^ tan ^ 

2 

i + e 


3 

[i + (e^-i)p(e)i 

2 e 



3 

(5 + 1)77(0-11 



e + 1 J 



(47) 

(48) 

(49) 


Herein, we assnme the simple case of a conformal ideal flnid, i.e. eiso = SPiso- As a resnlt, 
Eqs. (43)-(46) describe the eqnation of state of an anisotropic system of classical massless 
particles with vanishing chemical potential. 


B. Initial conditions 


In order to solve the set of partial differential eqnations (24)-(28), one has to specify 
the initial conditions at the initial longitndinal proper-time for the hydrodynamic evoln- 
tion, r = To, i.e. one has to define five three-dimensional prohles: A(ro, x_i_, <^), ^(tq, x_i_, <^), 
Ma;(ro, xj_, <^), Mj^(ro, xj_, <^), and 'd(ro, xj_, <^). During the initial moments of a heavy-ion colli¬ 
sion, due to inelastic interactions, the participating nucleons deposit energy into the space- 
time volume of the hreball. In this work, we assume that the distribution of deposited energy 
is well described by the optical Glauber model. As a result, the transverse momentum scale 
is given by 


A(ro,x_L,?) = £ 


-1 

iso 


^ 0 " 


p{b, x^,^) 


(50) 


p(0,0,0) 

where the proportionality constant eo is chosen in such a way as to reproduce the total 
number of charged particles measured in the experiment, and denotes the inverse function 
of eiso(A). For a central collision, in what follows, we identify Aq as the initial transverse 
momentum scale at the center of the simulated region. 

The density of sources is constructed using a standard mixed model 


p(5,x_l,<?)= (1 - +Pwn(^’X-l)) + 2kPbc(^,x_l) - ^ 5 ( 6 ,x_l)) , (51) 


where is the density of wounded nucleons from the left/right-moving nuclei and pbc 
is the density of binary collisions, both of which are obtained using the optical limit of the 
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Glauber model 




wwj 

1 —e A 2 J 


Pbc(&,x_l) = (TnnT + 




The longitudinal profile is taken to be [51] 

(<t - A<^)2 


h{q) = exp 


2a2 


e(|c|-Ac) 


(52) 

(53) 

(54) 


For the LHC case studied here, we use k = 0.145 for the mixing factor and an inelastic 
cross-section ctnn = 62 mb. The parameters of the longitudinal profile (54) were htted 
to reproduce the pseudorapidity distribution of charged particles with the results being 
= 2.5 and = 1.4. The shift in rapidity is calculated according to the formula [52] 


= 1 In + PwN + '^pjPwN Pwn) 

2 PwN + dwN ~ '^p(PwN ~ Fwn) 


(55) 


where all functions are understood to be evaluated at a particular value of b and x^. The 
participant velocity is dehned as vp = \/{^/s/2y — {niN /2Y/ {y/s/2) and rriM is the nucleon 
mass. In Eqs. (52)-(53) we have made use of the thickness function 


T{^p) = j dzpws(x±,z), 

where the nuclear density is given by the Woods-Saxon prohle 

— R 


(56) 


Pws{^±,z) = Po 


1 -|- exp 


-1 


(57) 


For Pb-Pb collisions, we use po = 0.17 fm ^ for the nuclear saturation density, R = 6.48 fm 
for the nuclear radius, and a = 0.535 fm for the surface diffuseness of the nucleus. 


V. RESULTS 

We now present our results for the photon spectrum and elliptic flow of QGP-generated 
photons. We focus here on Pb-Pb collisions with nucleon-nucleon center of mass ener¬ 
gies of y^NN = 2.76 TeV. Since the hard and soft contribution of differential photon rate 
described by Eqs. (8) and (9) are independent of the model assumed for the space-time 
evolution of the system, we hrst numerically compute the dimensionless differential pho¬ 
ton rate E/K^dR/<Rq on a uniformly-spaced three-dimensional grid in 0.1 < g^/A < 30, 
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47r?7/s^^^^^ 

0 

10 

100 

1 

0.552 

0.765 

1.009 

2 

0.546 

0.752 

0.992 

3 

0.544 

0.748 

0.990 


TABLE I. Values of the initial transverse momentum scale Aq in GeV used in all figures in the 
results section. 

0 < I 2 /I < 10, and —1 < log;^o('C + 1) < 2.5. The grid spacings used for these three variables 
were chosen in such a way that the full three-dimensional function is well approximated at 
continuous values using a spline-based interpolating function. For the hnal integration, we 
use the Vegas Monte-Carlo method to numerically integrate Eq. (17) over space-time and 
transverse momentum angle.® 

In all plots we assume a min-bias collision with b = 9.5 fm and we begin the aHydro 
evolution at tq = 0.3 fm/c. At tq, we assume that the produced matter has no transverse 
flow, i.e. Ux{To,^±,<i) = x_L, <^) = 0, while the initial longitudinal flow is of Bjorken 

form 'd(ro,xx,<^) = <^. We also assume that the initial anisotropy held is homogeneous, 
^(■^ 0 ) xx, <^) = ^ 0 - Finally, the initial central transverse momentum scale Aq used for all results 
is specihed in Table I. These values were tuned by requiring the hnal particle multiplicity 
to be a constant as t]/s and ^0 were varied. In all plots included herein we show results for 
central rapidity (?/ = 0) and have used a hxed = 0.32. 

We begin with Fig. 3, which shows the spectrum of medium photons obtained by integrat¬ 
ing over the full (3-|-l)D evolution of the QGP. The three diherent panels (a), (b), and (c) 
show the results obtained assuming initial anisotropies of .^0 = 0, 10, and 100, respectively. 
In each panel of Fig. 3, the lines correspond to three diherent values of the shear viscosity 
to entropy density ratio Aurj/s = 1, 2, and 3. As can be seen from panel (a), for .^0 = 0 
the resulting spectrum is nearly the same for all three values of 77/s, with only a very slight 
enhancement seen at large q±. The ehect of varying 77/s is larger in panels (b) and (c) in 
which the QGP was assumed to have an initially oblate momentum-space anisotropy. For 
^0 = 100, there is approximately a factor 2.5 variation in the medium photon spectrum at 
® During integration we set the rate to zero outside of the interpolated region. The excluded regions give a 
negligible contribution to the integrated rate. 
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FIG. 3. Medium photon spectrum for three different values of the initial anisotropy: (a) = 0, 

(b) = 10, and (c) .^o = 100. In each panel, the lines correspond to three different values for the 

shear viscosity to entropy density ratio dyrry/s = 1, 2, and 3. 

q± = Q GeV when varying rj/s between one and three times the lower bonnd. For q± ^ 2 
GeV, we see very little effect from varying rj/s for all valnes of considered. 

In Fig. 4 we present the same results in a slightly different manner. In this case, in 
panels (a), (b), and (c) we fix the shear viscosity to entropy density ratio to be Airri/s = 1, 
47 Tr]/s = 2, and Anri/s = 3, respectively. In each of the panels of Fig. 4 the lines correspond 
to three different values for the initial anisotropy in the system = 0,10, and 100. As can 
be seen from Fig. 4, there is significant variation in the high-energy photon spectrum as one 
changes the initial anisotropy of the QGP. For all values of rj/s considered, at g_L = 6 GeV, 
one finds that the QGP photon spectrum varies by approximately an order of magnitude 
when varying the initial anisotropy in the range shown. We note additionally that for qj_ < 


19 















(a) (b) 



(c) 



FIG. 4. Medium photon spectrum for three different values of viscosity (a) = 1 , (b) 

Aiirils = 2 and (c) = 3. In each panel, the lines correspond to three different values for the 

initial anisotropy in the system = 0,10, and 100. 

2 GeV, we see very little effect from varying the initial anisotropy. Taken together, we see 
that the low-energy photon spectrum is not sensitive to either or rj/s, if one keeps the 
hnal particle multiplicities hxed. 

In Fig. 5 we present the medium photon spectrum for six different values of the photon 
rapidity. As this hgure demonstrates for central rapidities, there is little dependence of pho¬ 
ton production on the photon rapidity; however, for > 3 there is a signihcant dependence 
on the rapidity. The fact that the central region is independent of the rapidity is consistent 
with the approximate boost-invariance of the quark-gluon plasma generated in heavy-ion 
collisions. The dependence on rapidity at forward/backward rapidity is due to both the 
breaking of boost-invariance in the realistic 3-|-ld anisotropic hydrodynamics code and also 
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FIG. 5. Medium photon spectrum for six different values of the photon rapidity. For this figure 
we assumed rj/s = 2/47r and = 0- 


the rapidity-dependence of the rate itself when the system is anisotropic (see Fig. 2). Note 
that, if the system was completely boost-invariant, then even with the rapidity dependence 
of the rate shown in Fig. 2, one would hnd that photon production does not depend on 
rapidity. We have verihed this explicitly. 

We now turn to the photon elliptic flow which results from the breaking of rotational 
symmetry around the beam-axis due to the development of hydrodynamic flow. To quantify 
the dependence of photon production on the momentum azimuthal angle, one makes a 
Fourier decomposition of the differential spectrum with respect to the momentum azimuthal 
angle 


dN 


dN 


1 -|- 2Vi cos (j)q + 2V2 COs(20q) -|- . . . 


(58) 


q±dq±dy d(j)q 27rq_i_dq±dy . 
where it is understood that the coefficients Vn are functions of q± and y.^ 

From the above relation we can compute V 2 in the usual manner by extracting the second 
Fourier coefficient from the series via 


V2{q±,y) = 




dN 


1 q^dq^dy d<pc 


COs(20q 


dN 


(59) 


q±dq^dy 

This coefficient is referred to as the ‘elliptic flow coefficient’, however, we emphasize that 
for photons a non-vanishing V 2 coefficient is not evidence of flow of the photons themselves, 
but instead the “imprint” of the transverse flow prohle of the QGP itself. 


® The coefficients are also implicit functions of the impact parameter, collision energy, colliding species, etc. 
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FIG. 6. Elliptic flow coefficient V 2 {q±, y = 0) for three different values of initial anisotropy parame¬ 
ter; (a) ^0 = 0) (b) '^0 = 10, and (c) = 100. In each panel, the lines correspond to three different 

values for the shear viscosity to entropy density ratio dvr/y/s = 1, 2, and 3. 

In Fig. 6, we show our results for the photon V 2 coefficient in three panels (a), (b), and 
(c) which show the results obtained assuming initial anisotropies of = 0, 10, and 100, 
respectively. In each panel, we once again show the result obtained for three different values 
of the shear viscosity to entropy density ratio of dvrry/s = 1, dvrry/s = 2, and Anrj/s = 3. 
As Fig. 6 demonstrates, we hnd that, regardless of the assumed initial momentum-space 
anisotropy, increasing the shear viscosity of the QGP results in an increase in photon V 2 - 
For the values of ry/s considered, we see at most a 300% increase in the photon V 2 with the 
maximum effect occurring at high transverse momenta. By comparing the three panels of 
Fig. 6 we also see that, for fixed t]/s, increasing the initial momentum-space anisotropy also 
results in an increase in the peak photon V 2 , however, as ,^o increases there is a reduction 
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FIG. 7. Variation of the photon V 2 with Tc. For this figure we took Airrijs = 1 and = 0- 

in the photon V 2 at large transverse momentum. In Fig. 7, we show the dependence of the 
photon V 2 on the assumed value for the transition temperature T^. In the hgure we show the 
results obtained assuming ATTr]/s = 1, .^o = 0, and G {175,155} MeV. As we see from this 
hgure, decreasing the critical temperature used when we integrate QGP emissions over the 
QGP four-volume results in an approximately 30% increase in the peak value of the photon 
V 2 - In addition, we see that the peak in V 2 moves to lower transverse momentum as Tc is 
decreased. 

Finally, as a cross check of the results shown thus far, in Fig. 8 we show (a) the photon 
spectrum and (b) the photon elliptic how as a function of the transverse momentum compar¬ 
ing the result obtained if we use the standard PMS criteria to set the soft/hard separation 
scale p* (red dashed line) or instead take the separation scale to be two times the PMS 
value (black line). As Fig. 8(a) demonstrates, using 2 p* results in a higher photon yield 
by approximately 30-50%. However, as Fig. 8(b) demonstrates, the additional production 
obtained when using 2 p* cancels in the ratio that determines V 2 and, as a result, the pho¬ 
ton V 2 is independent of the choice of the separation scale within numerical uncertainties. 
Finally, we mention that, although we only show the case = 0 and dvr/y/s = 3 in Fig. 8, 
for all ^0 and p/s considered here, we hnd that V 2 does not depend on the separation scale. 
The small differences between the two curves in Fig. 8 can be attributed to statistical errors inherent in 
the Monte-Carlo integration method used to evaluate the final spectrum. 
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FIG. 8. Dependence of (a) the medium photon spectrum and (b) the elliptic flow coefficient on 
transverse momentum for two different values of separation scale: p* and 2p*. In both panels, we 
assumed that the = 0 and 4:nr]ls = 3. 

VI. DISCUSSION, CONCLUSIONS, AND OUTLOOK 

In this paper we computed the photon spectrum and elliptic flow coefficient associated 
with photons emerging from the QGP as a function of transverse momentum. For the 
rate, we included the two leading-order processes necessary: Compton scattering and quark- 
antiquark annihilation. In order to properly deal with the IR divergences encountered in the 
rate calculation, we computed the soft and hard contributions separately. In the soft sector, 
the UV divergence was regulated by introducing a UV cutoff p* and, in the hard sector, 
the IR divergence was regulated by introducing an IR cutoff p*. As shown analytically by 
Braaten and Yuan [30], in the weak-coupling limit, the sum of the hard and soft contributions 
is hnite due to a cancellation between the UV and IR divergences and the result does not 
depend on the separation scale p*. However, for realistic couplings, one must evaluate the 
necessary integrals numerically. In this case, there is a residual dependence on the choice 
of the separation scale. To £x it, we used a PMS criteria to set p* to be the scale at which 
the derivative of the rate vanishes, thereby minimizing the dependence of the rate on the 
separation scale. 

After determining the rates, we then integrated them over the space-time volume of the 
QGP using (3-|-l)D anisotropic hydrodynamics to provide the space and time dependence 
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of the anisotropy ^(x), the transverse momentnm scale A(x), and the flow velocity u^{x). 
In the aHydro framework, the one-particle distribution function is guaranteed to be greater 
than or equal to zero unlike in the standard viscous hydrodynamic framework where there 
can be regions in phase space where the one-particle distribution function is negative. Our 
hnal results indicate that, if one holds the hnal particle multiplicity hxed, there is only a 
weak dependence of the photon spectrum on the assumed value of rj/s. However, we found 
that, for hxed r]/s, varying the initial momentum-space anisotropy .^o resulted in signihcant 
variations of the high-transverse-momentum photon yields. We found that, at = 6 
GeV, there was approximately an order of magnitude variation in the photon yield when 
varying 0 < .^o < 100. This offers some hope to constrain the degree of QGP momentum- 
space anisotropy by htting thermal plus prompt photon production at high energies to 
experimental data. We note that this is similar to the conclusion reached recently for 
dilepton production, where an enhancement of the high-energy dilepton production was 
observed when the initial condition was initially oblate in momentum space [53]. 

In addition to presenting results for the spectrum, in this paper we also calculated the 
elliptic flow coefficient V 2 associated with the azimuthal variation of the photon yields. Our 
results indicate that, for Attti/s = 1, the maximal photon V 2 coming from the QGP phase 
at LHG energies is approximately 2-3%. We hnd that increasing rj/s or results in an 
increase in peak photon V 2 . Our hnding that V 2 increases as rj/s is increased seems to be 
in conflict with some earlier papers, e.g. [8, 14-16], which found that incorporating viscous 
corrections resulted in a decrease in photon V 2 - One possible explanation for the different 
trends in the photon V 2 seen using aHydro versus second-order viscous hydrodynamics is 
that herein we employed the leading-order spheroidal form for the LRF distribution function, 
allowing for only one dissipative correction, quantihed by which maps to the longitudinal- 
transverse pressure anisotropy. In second-order viscous hydrodynamics, the description of 
the shear tensor is more complete, i.e. transverse pressure anisotropies and off-diagonal terms 
are present, resulting in a total of hve independent degrees of freedom. For quantitative 
assessment of photon V 2 , the additional viscous corrections could be important. Finally, 
we also note that our study is quite limited since we did not include fluctuating initial 
conditions or hadronic sources of photons. These are both major shortcomings of this 
work. Our intention herein was to study the systematics in a simply context in which the 
anisotropic screening and one-particle distribution functions were both taken into account 
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self-consistently. We plan to include multiple anisotropy parameters, fluctuating initial 
conditions, and hadronic emissions in forthcoming papers. Importantly, however, we point 
out that the dependence of the photon spectrum on the initial anisotropy found here is 
primarily sensitive to early-time longitudinal-transverse pressure anisotropies with the other 
viscous corrections being subleading. Therefore, we are confident that this effect is generic 
and reasonably well-described using leading-order aHydro. 

We mention in closing that, during the analysis, we showed that there is an approximately 
30% increase in the QGP photon spectrum when varying the hard/soft separation scale by 
a factor of two, however, we found that the photon V2 did not depend on the choice of the 
separation scale. In the future, in addition to improving upon the aHydro assumptions, we 
will combine the results obtained here with estimates for prompt photon production in order 
to extract constraints on early-time momentum-space anisotropies in the QGP. 
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